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Abstract. - We consider the surface critical behaviour of a semi-infinite two-dimensional 
layered Ising model, where the couplings perpendicular to the surface follow the aperiodic 
Rudin-Shapiro sequence. The model has unusual critical properties: depending on the strength 
of the modulation, at the bulk critical point the surface magnetization is either discontinuous 
or vanishes with an essential singularity. The critical surface magnetization as well as the form 
of the essential singularity are calculated exactly. 



Introduction. - Motivated by the discovery of quasi-crystals ||l| and the possibility to build 
artificial layered structures (see for example ) , the theoretical study of phase transitions and 
critical phenomena in quasi-periodic or, more generally, aperiodic systems has considerably 
developed in the last years (for a review, see [||). 

Following a series of numerical [4-9] and analytical [10-18] studies, an extension to aperiodic 
perturbations of the Harris criterion [|l^ for disordered systems was recently proposed [ po| . 
According to the Luck criterion, an aperiodic modulation can be irrelevant, marginal or 
relevant, depending on the value of a crossover exponent involving the correlation length 
exponent ly and the wandering exponent w, characteristic of the sequence. The form of the 
singularities for relevant modulations and the case of anisotropic critical systems were discussed 
in [^l| using scaling arguments. The criterion, first established for uniaxial aperiodicities, was 
extended to higher-dimensional ones in [|2^ . 

In this work we present exact results for the case of a relevant aperiodic modulation of 
the interactions following the Rudin-Shapiro (RS) sequence |2^, in the layered 2c? Ising 
model. Thus we complete a previous study [ ^ where exact results were obtained for irrelevant 
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and marginal perturbations, whereas the relevant case, which is the most difBcult to handle 
analytically, was only treated numerically and through scaling considerations. 

We consider the Ising model on a semi-infinite square lattice in which the couplings parallel 
to the free surface are constant and the modulation is associated with the perpendicular 
couplings. The problem is treated in the extreme anisotropic limit [ p5| where the transfer 
operator along the surface involves the Hamiltonion of a spin-^ Ising quantum chain 

n = --J2[Xk<y^,a^,^, + <yl]. (1) 

The surface magnetization can be obtained using standard methods p6| , as a function of 
the aperiodic couplings Afc with 

-1/2 

(2) 




Let us associate with each bond a digit fk which can take two values, or 1, depending on the 
type of the bond along the aperiodic sequence and let Xt = Xr^'' . The surface magnetization 
in eq. (0) can then be rewritten as 



-1/2 oo j 

, 5(a;,r)=^x^rJ-2"^ n^^^A, no = 0. (3) 

j=0 fc=l 



The critical coupHng Ac is such that limj^oo(l/i) X]fc=i -^fe — j2^, so that one obtains 
Ac = r~P°° , where poo is the density of 1 in the infinite sequence. 

Substitution matrix and scaling considerations. - An aperiodic sequence can be constructed 
through iterated substitions on the letters A, B, . . . . The RS sequence itself follows from 

A S{A) =AB, S{B) ^ AC, C -> S{C) ^ DB, D S{D) ^ DC. (4) 

Informations about the properties of the sequence are contained in its substitution matrix p9| 
with entries giving the number of letters i in S{j) {i,j = A,B, . . .). The length of the 

sequence after n iterations is related to the largest eigenvalue Ai through L ~ A" and the letter 
densities to the corresponding eigenvector. 

With the fc-th letter in the sequence, one associates a coupling which is either A or Ar, 
according to the value of fk (0 or 1). The cumulated deviation A{L) from the average coupling 
A at the length scale reached after n iterations can be expressed as 

L 

A{L) = Y,i^k - A) = A(r - l){nL - Lp^) ^ 5|A2|" ~ 6L^ , (5) 

k=l 

where 5 ^ r — \ and the wandering exponent uj |29| is given by ln|A2|/lnAi, where A2 
is the next-to-leading eigenvalue. Looking at the behaviour of the thermal perturbation 
A{L)/L ~ dL'^~^ under a change of scale by 6 = L/L', one obtains the scaling behaviour 
for the amplitude |l[ 

S' = b'^/'^S, <i> = l + u{uj- 1). (6) 

The relevance of the perturbation depends on the sign of the crossover exponent For the 
2d Ising model with 1^ = 1, the aperiodic modulation is irrelevant when w < 0, marginal when 
uj = and relevant when a;>0. In the marginal case, varying exponents are obtained p^ . 



F. iGLOi et al: relevant aperiodic modulation 



93 



The behaviour of the system in the vicinity of the critical point follows from scaling 
considerations. The surface magnetization behaves as 

ms(i,(5) =<AF(u), u^j, (7) 

where t = 1 — (Ac/A)^, /3s is the magnetization exponent at the ordinary surface transition, 
^ ~ t~'^ is the bulk correlation length and the scaling variable u involves a new characteristic 
length associated with the aperiodicity, Za ^ IJI^"^/*, which stays finite at the bulk critical 
point. 



Surface magnetization for the Rudin-Shapiro sequence. - For the RS sequence, one may 
associate fk=0 with A and B and fk = 1 with C and I? in (^, starting the sequence on A. 
The largest eigenvalues in modulus of the substitution matrix are Ai = 2, A2=±V2, leading 
to the wandering exponent 'jJ = \- The asymptotic densities of the letters yield p^o — \ thus 
the critical coupling is Ac = r^^^"^. 

Using the two-digit correspondence A, B, C, 13 — > 00, 01, 10, 11 which follows from after 
one iteration, one obtains the following recursion relations for the n^'s: 



n-8p+2 = 2p + 2n2p+i, 
n-8p+4 = 2p + 1 + 2n2p+i, 
n-8p+6 = 2p + 1 + 2ri2p+2, 

This allows us to rewrite the sum in (||) as 

S{x, r) — —x'^ 



mp+i = 2p + 2n2p + /2p+i, 

"8p+3 = 2p + 2n2p + 3/2p+i, 

«8p+5 = 2p + 1 + 2n2p+l + /2p+2, 

«8p+7 = 2p + 2 + 2n2p+l + /2p+2- 



(8) 



1 + :e + -;-{r + r + 



2x4 



r + r 



r 

„2\ 



S{x\r^) 



(9) 



where x — (Ac/A)^ 

r^). At the critical point x 
functional equations: 



Changing r into — r leads to a second equation giving 5(x^,r^) and 
1, with S{l,r) = l + (7(r), one obtains the following 



1 



9{r') ^ -A9{r) + g{-r)] , g{-r') 



(1 -r)4g(r) + (1 + r)4g(-r) 



4(l + r-2)2 



(10) 



Now, when r < 1 one may look for g{r) under the form of a series expansion in powers of r, 
g{r) = ^^Qdpr^, where, according to (|l^), the coefficients satisfy 



a2p — 2ap, 



-1 + a4p_3 — 4a2p-i, a4p+i + a^p^i — a2p+i 



These recursion relations are solved by Op =p, leading to g{r)~r/{l~r 
to (S), the critical surface magnetization is given by: 



1 



(r<l). 



2a2p + a2p-i. (11) 
l'^ and finally, according 

(12) 



When r> 1, g(r) diverges, i.e., the surface magnetization vanishes at the bulk critical point. 
To study its behaviour on a finite system with size L, one has to terminate S{l,r) in (^ at 
j — L. For L = 2^"~^ one finds by inspection that it corresponds to the first 2" terms of the 
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power series for g(r), up to a correction for the last one. Thus, one obtains 



"-1^2" 



r - 1 



1 



2r 



-2" - 1 



(13) 



When n ^ oo the previous result is recovered for r < 1, whereas the expression diverges for 
r>\. In the scaling limit, r— l<Cl, L—fOo, one obtains the following finite-size behaviour for 



the critical surface magnetization: 



"r- 1' 


1/2 


.V2L. 


exp 



(r - 1)V2L 



(14) 



Thus, the finite-size dependence of TOs,c, having a stretched exponential form, is anomalous 
for r>l. A similar behaviour is obtained below for the temperature dependence. 

In order to determine the general behaviour of the surface magnetization when 1 ^ t > 
and |r— 1| <C 1, we take a continuum approximation for the sum in (|^) and use (^ to write 

j-n^jl Poc-Aj'^ , so that 



/ AjeM-jt + Bf): B = A 
Jo 



Inr ~ A{r - 1). 



When B <0, expanding exp{—jt) and integrating term by term, one obtains 



■E- 

Tl = 



n + 1 



r(i/^) 



r(2/a>) 



r(l/a;)|B|i/- 



(15) 



(16) 



The surface magnetization approaches its nonvanishing critical value as ms{t)—ms^c^t^'' , with 
a critical exponent /Jg = 1 for any value of the wandering exponent uj. For the RS sequence 



with (^=2, one obtains 



Ail 



V2 



3t 



A2(l-r)2 



(17) 



in agreement with the exact result in (|T^) and with the scaling form in (0) with /3s = i , ly—l 
for the 2d Ising model, and i^(u) ~u~^/^(l + au + ■ ■ ■). According to (|^), the critical surface 
magnetization generally vanishes at r = l as |r— 1|'^=/*. 

When B>0, the integral in ( |l5| ) can be evaluated using Laplace's method when 0<a;<l. 
A straightforward calculation then gives, up to a power law prefactor. 



exp 



[uj - 1) 



2w(tjB)i/('^-i) 
where the last expression is the RS result. 



exp 



A^ (r - 1)2 



(18) 



Discussion. - The surface order ( |12D at the bulk critical point is linked to a local 
enhancement of the couplings for r< 1. Starting the sequence on D would amount to exchange 
weak and strong couplings and the surface first-order transition would then occur for r > 1 . 
Since Poo = ^, one does not expect such an asymmetry in the bulk critical behaviour. 

The amplitude A in ( p^ , which was considered as a constant in the continuum calculation, is 
actually an oscillating function of j ||2^ . This does not affect the stretched exponential in ( p^ ) 
but may influence the prefactor which was ignored. In order to determine it, one may use the 
exact finite-size result (ph, replacing L by the relevant thermal length L, associated with the 
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aperiodicity, which is obtained by equating the two terms in the exponential of eq. ([15|). Since 
[(r— — [(r—lj/t]'^ for the RS sequence, the prefactor is proportional to t^^^ and 
the scaling function in eq. ([7|) is simply F{u)^exp{—bu~^). 

To conclude, we would like to emphasize the close relationship between our results and 
those obtained for relevant extended perturbations |2^, In these systems, the couplings 
deviate from the bulk ones by AL~y , where L is the distance to the surface. For the aperiodic 
systems, this corresponds to the average deviation per bond A{L)/L^L'^^^ at a length scale 
L. Relevant perturbations are obtained for y <1, which corresponds to w > 0. Close to the 
critical point, the short-distance structure of the perturbation is irrelevant and the surface 
singularities are the same for both systems. 
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